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Abstract 

Neutrino oscillation experiments often employ two identical detectors to 
minimize errors due to inadequately known neutrino beam. We examine var- 
ious systematics effects related to the prediction of the neutrino spectrum in 
the 'far' detector on the basis of the spectrum observed at the 'near' detector. 
We propose a novel method of the derivation of the far detector spectrum. 
This method is less sensitive to the details of the understanding of the neutrino 
beam line and the hadron production spectra than the usually used 'double 
ratio' method thus allowing to reduce the systematic errors. 
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1 Introduction 



Neutrino oscillations experiments can be divided into two classes: 

• disappearance experiments 

In these experiments one measures a deficit of the observed rate of neutrino 
interactions with respect to the rate expected from the knowledge of the neu- 
trino beam. 

• appearance experiments 

In these experiments one detects interactions of the neutrino flavor not present 
in the neutrino beam at its source 

A positive result of an appearance experiment constitutes an undisputable proof 
of the existence of the neutrino oscillations and it is fairly independent of the detailed 
understanding of the neutrino beam. This is often contrasted with the disappear- 
ance experiments, where the uncertainty of the neutrino flux at the position of the 
detector usually dominates the systematic error. 

Further steps of the determination of the oscillation parameters involve a quan- 
titative estimate of the rate of appearance or disappearance and their dependence 
on the neutrino energy. At this stage the understanding of the neutrino flux at the 
position of the detector constitutes a major source of the systematic error. 

Absolute prediction of the neutrino fluxes is difficult, as the saga of measuring of 
the total cross sections teaches us. Modern tools of detailed simulations make the 
predictions somewhat more reliable, but accuracy of these predictions is limited by 
the scarcity of the data constraining the shower simulation code. 

A powerful method of improving the knowledge of the neutrino flux at the posi- 
tion of the detector consists using two detectors in the same beam line: 

• a 'near' detector positioned near the neutrino source to establish the beam 
characteristics and to enable a reliable prediction of the neutrino flux at the 
'far' detector 

• a 'far' detector positioned at the distance where oscillation effects are expected 
to be present. 

This technique was pioneered by the CCFR|I| collaboration at Fermilab and then 
subsequently used by CDHSWg and CHARM § experiments at CERN. Recently 
the two detector approach was employed by the K2K j|] experiment in Japan. 

High statistics of the neutrino interactions expected at MINOS offer a potential 
for precise determination of the oscillations parameters. This, in turn, requires a 
very good understanding of the neutrino flux. 

This note describes contributions to the systematic errors of the neutrino flux 
prediction and proposes a novel method of predicting the neutrino flux with reduced 
sensitivity to various systematic effects. 
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2 Neutrino Beam 



Neutrino beam is produced in three steps: 

• pions and kaons are produced in the target 

• (optional, but necessary to attain high beam intensities) produced particles 
are focused directed towards the detector 

• pions and kaons decay inside the decay volume producing neutrinos 

We examine these steps at some level of details to indicate the possible sources 
of systematic errors of the neutrino flux predictions. 

2.1 Particle production in the target 

Primary protons strike the target and produce pions and kaons. For low energy 
beams and long baseline experiments (hence neutrinos produced at small angles) 
pions produce the main contribution to the neutrino flux. 

Production spectra of pions produced by protons are known rather poorly. Sur- 
prisingly large uncertainties exist in the knowledge of the inclusive cross sections, 
even for the elementary pp reaction. Situation is further complicated by the 
fact that the neutrino beams utilize extended targets, where the absorption and 
re-interaction of the particles produced in the first interaction play a significant 
role. 

Uncertainty in the production spectra is a direct consequence of a scarcity of the 
relevant experimental data. The primary source of the experimental information 
is the measurement of the particles yields in pBe collisions at 400 and 450 GeV/c 
performed at CERN as a tool for understanding the CERN NBB and WBB neu- 
trino beams', ||], |7j. Extrapolation to different target materials is guided by the 
measurement of Barton et al.|| 

There are two approaches to predictions of the particle fluxes for the neutrino 
beam: 

• use the existing data to constrain and tune complete shower simulation codes. 

This approach is represented by FLUKA |{J and MARS [[Kj. It should be 
pointed out that we use FLUKA as implemented in GEANT (often referred 
as GFLUKA) . There exists a more modern version of FLUKA, re-tuned using 
the SPY|| [7| data, but unfortunately it is not interfaced to GEANT yet. 

• parameterize the existing data with an ad-hoc analytical formula and use it as 
a particles source function. Some physics-driven scaling laws must be employed 
to enable extrapolation of the data to different targets, proton energies and 
phase space regions of the produced particles. 
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This approach is represented by MalensekjT]]] and BMPT[|TI| 'models'. The 
BMPT parameterization covers a complete phase space of the produced par- 
ticles. It is based on much more complete data set and it is developed for 
predictions of neutrino beam produced by different energy protons. Malensek 
formula is an interpolation of the Atherton|| data derived for a specific case of 
the Fermilab NBB. Use of this formula outside the region covered by the data 
points used is ill justified. We use it here, anyway, as an example of rather 
extreme variation of the production spectra. 

The primary difference between various production models is the absolute flux 
predictions. A typical range of the variation is of the order of 20%. Whereas this is 
an important practical factor, it is not important for the prediction of the shape of 
the spectrum observed at the far detector. 

A potential source of the systematic error on the shape of the far detector spec- 
trum is related to the shape of the p t distribution of the produced particles. The 
shape of the p t determines the angular divergence of the hadron beam after the fo- 
cusing elements and consequently the distribution of the decay points of pions along 
the decay volume. 
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Figure 1: Invariant cross section for ix + production as a function of px- Lines 
represent prediction of the BMPT model for different longitudinal momenta of 7r + . 
Data points are from the Refs. [^|, |5]. [7|] 



Figs.[TJ and ^ compare the shape of the p t spectra predicted by the BMPT and 
MARS models |13]]. They show that a considerable uncertainty exists in the shape: 
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physics- inspired general models tend to have an exponential fall-off with p t , whereas 
the experimental parameterization follows the experimental data better. 
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Figure 2: Invariant cross section for 7r + production as a function of px- Lines 
represent prediction of the MARS model for different longitudinal momenta of 7r + . 
Data points are from the Refs. |5], |6], [7[] 

A relationship between the p t distribution of produced pions and the neutrino 
spectra is illustrated in Fig.|] by showing the distributions of p t of pions responsible 
for neutrinos of different energies, as predicted by different production models. 

Systematic errors due to hadron production rates on the predicted neutrino fluxes 
have been extensively studied recently in general context JT3| , as well as in the context 
of NuMI/MINOS experiment [0, 

2.2 Focusing Elements 

A system of parabolic magnetic horns focuses the secondary particles by giving them 
a p t kick proportional to the radius at which particle traverses the horn. In case 
of a pointlike source such a system focuses all particles produced with a particular 
momentum. Finite radial sizes of the horns lead to a finite momentum byte of 
the secondary particles which are focused. The central momentum of the focused 
particles depends on the relative distance of the target and the horn. 

The hadron beam does have a significant angular divergence due to several ef- 
fects: 
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Figure 3: Distributions of the production p t of neutrino parents weighted by the 
neutrino flux detected in the near (left column) and far (center column) detectors, 
predicted by different hadron production models. Ratio of the p t distributions pro- 
ducing neutrinos of different energies in the far and near detectors is shown in the 
right column. 

• particles produced at small angles pass through the opening of the magnetic 
horns and enter the decay volume with the angular divergence characteristic 
for the bare target beam. This component is defined by the aperture of the 
horn opening and the distance of the target from the horn. These particles 
produce a high energy tail of the neutrino beam. This component produces 
also the main contribution of the background. 

• finite size of the primary proton beam 

• the finite size of the production target. Particles produced upstream of the 
focal point the the first horn will be overfocused whereas particles produced 
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downstream of the focus will be underfocused, thus producing a divergent 
beam. 

An overall angular distribution of hadrons after the focusing system is shown in 
the Fig.|. 
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Figure 4: Angular distribution of neutrino parents at the entrance to the decay pipe. 
2.3 Finite target length 

The target length is typically comparable to the distance between the target and the 
first horn. As the result, different parts of the parent momentum spectrum will fall 
into acceptance of the horn for different sections of the target: higher momentum 
particles from the the beginning of the target and lower momentum particles from 
the end of the target will be focused. This will result in the difference in the spectra 
of the resulting neutrinos, as shown in the Fig.[5| The three effects are readily 
apparent: 

• A number of pions (and hence a number of neutrinos) produced decreases 
along the target. This is due to the attenuation of the primary proton along 
the target. 

• Spectrum of pions (and neutrinos) is shifted towards lower energies along the 
target. This is the effect of the finite acceptance of the horn focusing system. 

• Relative contribution of the high energy tail increases along the target. This 
tail is due to the small angle particles which pass through the opening in the 
horn. The angular acceptance is much smaller for particles produced in the 
front section of the target. The attenuation of small angle particles traversing 
the target contributes another suppression factor. 
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Figure 5: Neutrino spectra from decays of particles produced in different sections of 
the target. 

The long production target in conjunction with the two- horn focusing system 
will produce the pion/kaon (neutrino) beam with characteristics dependent upon the 
hadron energy. Fig.|6| shows the distribution of the angle of the neutrino parents after 
the focusing system. Low energy part of the spectrum comes from well focused beam. 
Neutrinos in the range 5 to 7 Gev are produced by pions which are increasingly less 
and less focused, whereas the high energy tail E u > 8 GeV are produced from pions 
produced at very low angles and passing through the opening of the magnetic horns. 

2.4 Decay Volume 

Pion and kaon beam formed by the horns decays in the decay volume consisting 
of the last section of the target cave and the decay pipe. Neutrinos are produced 
in two-body decay n + — > n + + z/ M . Energy and the flux of the produced neutrino 
depend on the decay angle 6 dec as: 

f 0A3E * m 



Flux 



27 



A 



(2) 
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Angle of the neutrino parent decaying in the decay volume 
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Figure 6: Pion angle at the decay point for different pion (neutrino) energies observed 
in the near detector. 



where 

• 7 = is the Loretnz boost factor of a pion 

• 9 is a decay angle, i.e. the angle between the pion and the produced neutrino 
directions 

• A is the area of the detector and z is the distance between the decay point 
and the detector. 

Finite transverse size of the decay volume, while important to maximize the 
overall neutrino flux, causes a systematic difference between the neutrino spectra 
observed at the near and far detectors. This is due to the fact that the decay angle 
(i.e. the angle between the parent pion and resulting neutrino directions) necessary 
to reach the near and the far detectors are different for decays occurring at finite 
radii as shown in Fig.[7[ (A necessary element of a rigorous proof is the approximate 
azimuthal symmetry of the hadron beam and the fact that ^ > for the majority 
of the focused particles.) 

Fig.|| illustrates the dependence of the produced neutrinos and Fig.|] shows the 
neutrino flux as a function of the decay angle. This dependence is very strong for 
high energy pions, say E n > 20GeV. At the energies relevant to the main component 
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Figure 7: Differences in the decay angles to reach the near and the far detectors due 
to finite decay pipe dimensions. A systematic difference A@ = ©at — ©f is increases 
with the z of the decay point. 



of the low energy beam, E n < 10 GeV neutrino energy is very weakly dependent on 
the decay angle for 9 dec < 2 mrad. 

Hadron beam produced by the focusing system has an angular divergence of the 
order of few miliradians. Such a divergence produces a hadron beam much larger 
than the radius of the decay volume, therefore most of the particles will hit the walls 
of the decay pipe, unless they decay earlier. Dependence of the angular distribution 
of pions on the pion energy (Fig. |6|) will therefore lead to a correlation between 
the resulting neutrino energy and the average decay point of the parent pions as 
shown in Fig.|TD|. This effect is contributing to a difference between neutrino spectra 
observed with the near and far detectors. The neutrino flux observed by the detector 
varies with as Q>(z) ~ -, — — where Zd ec is the average position of the decaying 
pion. This reduction of the flux due to the angular divergence is independent of the 
neutrino energy, as the variation of Zd ec is negligible compared to [zdet — Zdec)- In 
case of the near detector, though, the variation of of the attenuation of the flux with 
the neutrino energy is significant. 

Finite transverse size of the decay volume, in conjunction with the angular diver- 
gence of the pion beam leads to another contribution between the neutrino spectra 
observed in the near and far detectors. For a pion decaying at some radius R inside 
the decay volume the decay angle pointing to the near detector is larger than the 
decay angle necessary for neutrino to reach the far detector, see Fig.|TT|. 

As a consequence of the dependence of the neutrino energy on the decay angle 
(FigH) the neutrino energies detected at the near detector will by systematicly lower 
in the near detector even if both neutrino fluxes are produced by the same parent 
hadron beam, Fig .[12. 
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Figure 8: Neutrino energy from the ir + — > /x + + decay as a function of the decay 
angle for different pion energies. 



2.5 Near and far neutrino fluxes 

Hadron beam produced at the target and focused by the horns gives rise to the 
neutrino beam. Neutrino fluxes detected by the near and far detectors are highly 
correlated, as they are produced by the same hadron beam. They are not identical 
though. The difference of the spectra at the far and near detectors are primarily 
due to two effects: 

• a solid angle subtended by the unit area at the near detector varies considerably 
between the beginning and the end of the decay volume. This effect is very 



small for the near detector. This effect is shown in Figs.[13] and [14] 

As the result, the neutrinos produced towards the end of the decay pipe con- 
stitute much larger fraction of the neutrino flux detected at the near detector 
in comparison with the far detector. As the distribution of the decay points 
along the decay pipe varies with the hadron momentum, the resulting energy 
distributions at both detectors will be somewhat different, leading to a varia- 
tion of the ratio Rf/n = -j# as a function of neutrino energy Thus the 
prediction of the far detector spectrum from the spectrum measured at the 
near detector relies on the proper modeling of the longitudinal distribution 
of the decay points. The shape of longitudinal distribution of decay points 
is primarily determined by the width of the decay volume and the angular 
divergence of the hadron beam. At low hadron momenta a finite lifetime of 
pions and kaons has a significant contribution. 
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Figure 9: Neutrino flux from the 7r + — > /i + + z/ M decay as a function of the decay 
angle for different pion energies. 



• angular divergence of the hadron beam spreads particles throughout the vol- 
ume of the decay pipe. Neutrinos from decays at large radii will be produced 
at different decay angles to reach the near and the far detectors. As the result, 
the energy of neutrinos observed at these two detectors will differ, with the 
energy at the far detector being in general higher, as shown in FigJT5|. 



3 A method to derive the Far Detector spectrum 

A neutrino spectrum observed at the far(near) detector can be derived from the 
knowledge of the secondary hadron distribution (after the focusing elements) and 
the geometry of the decay volume: 



dN v 

_ near 

~d~E 



and 




dE 



far 



where: 




Fk/k (ri,9, p) P n /K (ri,8,p,z) W n/K ( 

(3) 

Fir/Kin, 9,p)P 7r / K (r i , 9,p, z)W 7r/K (z, r dec , 6,p, z far , E far ) dr { dOdpdz 

(4) 
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Figure 10: Decay point of the neutrino parent for different neutrino energies. 



• F n /K(ri,9,p) is the radial, angular (6 = ^Jp\ + Vy/Vz) and momentum distri- 
bution of pions(kaons) after the focusing elements. This function is a convo- 
lution of the production cross section, horn acceptance and horn focusing. 

• P-7r/K(ri,9,p, z) is a probability that a pion(kaon) with momentum p, radial 
position Ti and the angle 9 will decay at the position z along the decay volume. 
The radial position of the decay will be r dec = + 9z. For simplicity we have 
set z = at the end of the focusing system. This function depends solely on 
the geometry of the decay volume 

• W n / K (z,r dec ,p,9,z near/far ,E near / far ) is a probability that a pion(kaon) with 
momentum p and the angle 9 decaying at the position z along the decay 
volume at the radius r dec will produce neutrino with energy E near /f ar at the 
center of the respective detector. This function is determined purely by two- 
body kinematics. 
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Figure 11: Decay angle to reach the near detector vs the one to reach the far 
detector. Spread of the angles is determined primarily the the transverse size of the 
decay volume. 



For simplicity of the presentation we have assumed an axial symmetry of the 
problem. In practice the integral is over a six-dimensional space of (x, y, p x , p y , p z , z). 
A complete treatment must also include neutrinos produced before the end of the 
focusing system. 

3.1 First Approximation: Pencil-like Beam 

In case of a pencil-like hadron beam, i.e. a beam with no angular divergence, 
directed towards the far and near detectors the relation between the near and far 
neutrino spectra is simplified. The decay angle 9 = and the neutrinos detected at 
both detectors have the same energy E near = Ef ar = 0A3E W . Thus the measured 
neutrino flux at the near detector provides a direct measurement of the energy 
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Figure 12: Energy of the neutrino observed in the far detector vs the energy observed 
in the near detector. 

distribution of the parent pion beam. 

The neutrino spectra observed at the near and far detectors are not identical, 
though. They differ, because of the differences in the acceptances of the near and 
far detector. The acceptance is a function of the distance from the decay point to 
the detector. This distance is practically the same for all neutrino energies in the 
case of the far detector, but it is a function of energy, due to the finite pion life-time, 
at the near detector position. Thus the neutrino spectrum at the far detector can 
be derived from the spectrum observed at the near detector as: 



dN 



far 



dE 



dN 



far 



dE 



dN n 



dE 



dN 



MC 



dE 



T{E) 



dNr, 



dE 



(5) 



The transfer function T(E), in this simple case, is in fact calculable from the 
known pion life-time, the length of the decay region and the position of the near 
and far detectors. 



18 




20000 40000 60000 
Z decay, 2enu3 gev 



x 10 



20000 40000 60000 
Z decay, 2enu3 gev 




20000 40000 60000 
Z decay, 3enu4 gev 




20000 40000 60000 
Z decay, 4enu5 gev 



x 10 

0.8 
0.6 
0.4 
0.2 




20000 40000 60000 
Z decay, 3enu4 gev 





20000 40000 60000 
Z decay, 4enu5 gev 




20000 40000 60000 
Z decay, 5enu6 gev 



20000 40000 60000 
Z decay, 5enu6 gev 



Figure 13: Decay point of the neutrino parent for different neutrino energies. On the 
left are decays yielding neutrino in the near detector, on the right there are decays 
yielding neutrino in the far detector 



3.2 Second Approximation: Small Aperture Beam and 'Dou- 
ble Ratio' Method 

Focusing of particles with large momentum spread and produced in an extended 
target into a parallel beam is not possible. A realistic beam has a very sizeable 
angular divergence ( see Fig||) which is different for pions of different energies. This 
effect determines the effective distribution of the decay points along the length of 
the decay volume, as shown in Figs.0 and 14. Derivation of the transfer function 



T(E) requires, in this case, a Monte Carlo simulation of the production and focusing 
of the pions and kaons. The transfer function T(E) is just a ratio of integrals [| and 
[3] and its shape is shown in Fig.|16| The mean value of the function, ~ 0.15 x 10~ 5 , is 
given by the square of ratio of the far and near detector distances from the neutrino 
source. The shape of this function is related to the variation of the longitudinal 



distribution of the decay points along the beam line as shown in Figs. 13 and 14 
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Figure 14: Decay point of the neutrino parent for different neutrino energies. On the 
left are decays yielding neutrino in the near detector, on the right there are decays 
yielding neutrino in the far detector 



The structures observed in this function are reflections of the acceptance of the 
focusing elements. A sharp rise around 5 GeV is a related to the fact that the pions 
with E n > 12 GeV start missing the first horn and loose benefit of the focusing. 
As a result their angular distribution broadens, and they hit walls of the decay 
volume before they decay. The average distance from the decay point to the near 
detector increases (in comparison with the focused pions of lower energy) and the 
flux observed at the near detector is reduced thus producing a rise in the Far /Near 
ratio. An analogous mechanism, involving the horn 2 is responsible for the rise of 
the Far /Near ratio around 8 GeV. 

The transfer function will depend, somewhat, on the assumed production spectra 
in the target, but its overall value and the detailed shape is chiefly determined by 
the geometry of the beam line. 

For a beam of a very small radial aperture, such that decay angles leading to 
the near and far detectors are very similar, the energy of neutrinos observed at the 
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Figure 15: Ratio of the energy observed in the far and near detectors for neutrinos 
produced at different radii in the decay volume 



both detectors will be very close (see Fig.|l5D and the dependence on the production 



dN 



far 



I dN near \ 
I dE ) j 



cross sections will cancel, to a first order, in a ratio , ,,. 

V dE ■ «"=- / MC 

This method of predicting the far detector spectrum, often referred to as a 'double 
ratio' method is widely used to estimate the systematic errors in the disappearance 
experiments. Owing to a smooth behavior of the underlying production cross sec- 
tions, acceptances and focusing, this is a quite robust method of predicting the far 
detector flux. 

Shortcoming of this method comes to light when some instrumental effects mod- 
ify the neutrino (or pion) spectra in a very limited energy range. An example of such 
an effect is a small displacement of focusing elements (horns) whereby some small 
part of the phase space of the produced particles becomes focused and produces neu- 
trinos of some energy, at the expense of particles in a neighboring part of the phase 
space, being moved out of the horns acceptance and hence producing significantly 
fewer neutrinos. As a result of the extended transverse size of the decay volume 
these 'gained' and 'lost' neutrinos will have somewhat different energies in the near 
and far detectors, thus leading to a characteristic bipolar shape of the predicted far 
detector spectrum fluctuating around the real observed spectrum. 

The limitation of the 'double ratio' method becomes apparent when one attempts 
to predict a spectrum of the far detector in the 'off-axis' position, where the spectrum 
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Figure 16: Ratio of the neutrino fluxes observed in the far and near detectors as 
a function the neutrino energy. Top figure is for all neutrinos, bottom figure is for 
neutrinos originating in the first section of the decay volume, Zdecay < 400m 

of the observed neutrinos is completely different then that in the near detector. 

3.3 A Far-to-Near Correlation Matrix Method 

Neutrino spectra in the far and near detectors are correlated, as they result from 
the decays of the same parent hadron beam. Thus, every neutrino with energy 
E n observed at the near detector implies a certain flux of neutrinos with energy 
distribution Ef given by 

dN = W(z,r dec ,p, z f ,Ef) 
dE f W(z,r dec ,p,z n ,E n ) 

This flux, implied for the far detector, depends on the position of the decay point as 
well as on the momentum vector of the parent hadron. Integration over all decays 
yields a matrix, which correlates the spectra at both locations. The matrix element 
M AEfAEn is given by 

/ • • • / F w (n, 9, p)P(r h 9, p, z)W(z, r dec , 9, p, z f , E f ) dr { d9 dp dz dE f 



MAE f AE n 



J ■■■ J F n {r h 9, p)P(r h 9, p, z)W(z, r dec , 9, p, z n , E n ) dr { d9 dp dz, dE n 



(7) 
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where the integrals are over the bin sizes AEf and DeltaE n , respectively. In 
order to relate the observed event spectra (rather than the neutrino spectra) at 
both detectors the weight functions W(z, rd ec , 6,p, z n /f, E n /f) should be replaced in 
the integral (|7|) by W(z, r dec , 9,p, z n/f , E n/f ) x a^ t (E n/f ). 

Given a vector containing the observed event spectrum in the near detector, 
iV = (Ni, N 2 , ■ ■ ■ , N n ) the predicted far detector spectrum F = F 2 , ■ ■ ■ , F n ) can 
be derived as: 









M 12 . 


■ M ln ~ 






(F 1 ,F 2 ,-- 


■,F n ) = 


M 2 i 


M 22 . 


■ M 2n 




N 2 






M nl 


M n2 . 


• M nn 







(8) 



We have computed the correlation matrix M for the event spectra in 0.5 GeV 
bins using a set of Ntuples generated using GNUMI simulation program and corre- 
sponding to 1.7 x 10 7 protons on target Uppermost 15 x 15 fragment of the 



matrix, thus showing the correlation between far and near detector event spectra 
below 7.5 GeV is shown in (||). Columns of the matrix are related to the energy 
observed in the far detector, rows to the energy in the near detector. 
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x 1(T 6 (9) 



The matrix provides a very good representation how the far detector spectrum 
relates to the near one. Column 10, for example, states that 1 event observed in 
the near detector in the energy bin between 4.5 GeV and 5 GeV implies (modulo an 
overall factor 10 -6 ) that the far detector should register 1 event with energy between 
4.5 GeV and 5 GeV, 0.47 events between 5.0 GeV and 5.5 GeV, 0.07 events between 
5.5 GeV and 6.0 GeV etc. 

The matrix is computed for all the simulated v a CC events, regardless of the 
actual neutrino parent. To a good approximation, no background from interactions 
of other neutrino types (including v^) needs to be included, as the detection of a /i~ 
uniquely identifies a v u charged current event. For the present studies, smearing was 
ignored, i.e., matrix elements were calculated from generated near and far detector 
energies. 

For the analysis purposes the correlation matrix M should be computed for the 
spectra with smearing and acceptance cuts properly included. 
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4 Systematic Errors: Sensitivity of the Prediction 
to the Details of the Parent Hadron Beam 



When the near-to-far correlation matrix M is computed in the same model as the 
near and the far detector spectra, the predicted far spectrum using the equation |8| 
matches exactly the actual far spectrum. It implies that the far detector spectrum 
would be known exactly, should the beam simulation program describe the real 
beam in all details. In practice, the understanding of the beam is somewhat limited, 
hence leading to systematic errors on the far spectrum prediction. In general, these 
errors will be due to the differences of the spatial and momentum distribution of 
pions decaying within the decay volume. 



4.1 How Small the Systematic Errors Should be? 

Clearly, smaller the systematic errors are - the better. In practice there is always a 
threshold below which a further reduction of the systematic error of the experiment 
has no impact on the physics results. 

A judgment of the impact of the systematic errors on the physics result is fairly 
complicated. It depends on several aspects of the experiment, like: 

• how large the measured effects are with respect to the systematic errors? 

• what is a type of a measurement? For example a fairly small, but systematic 
distortion of a spectrum may have a significant influence on the determined 
shape parameters. 

• what are the experimental smearing effects? 

• what are the statistical errors? etc... 

Neutrino disappearance experiment, like MINOS, is expected to detect and mea- 
sure a very dramatic effect of 100% disappearance at a particular neutrino energy. It 
is very difficult to imagine systematic errors related to the understanding of the neu- 
trino beam which could create false effects of this magnitude or could significantly 
distort the observed effect of full disappearance. 

On the other hand a detailed shape of a disappearance curve has a significant 
physics interest. In the least interesting case it will provide yet another proof of the 
underlying physics mechanism. In a more interesting case it may provide a clue for 
physics phenomena beyond the current orthodox picture. 

Statistical errors of the experiment provide a reasonable yardstick for the sys- 
tematics error. The expected statistical errors of the observed neutrino flux at the 
5 kton far detector and an exposure corresponding to 8 x 10 20 protons on NuMI 
target are shown in Fig. 17 for the case of 0.5 GeV bins. These error refer to the 
'no-oscillation' case. 
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Statistical errors, 0,5 GeV bins 
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Figure 17: Statistical errors in 0.5 GeV bins of far detector neutrino event spectrum 
for 5 kton detector and exposure equivalent to 8 x 10 20 protons on target 



In reality the reduction of the statistics induced by the neutrino disappearance as 
well as the experimental smearing effects will tend to further reduce the importance 
of the systematic errors. 

4.2 Dependence on the Particle Production Model 

The extent, to which the distributions of particles in the integral (^) defining the 
correlation matrix change with different production models, and, as a result, the 
matrix elements change, contributes to the systematic error. 

Dependence of the predicted far detector flux can be evaluated by applying 
the same correlation matrix M to the near detector spectra, as derived in various 
production models, and comparing the predicted far detector spectra computed from 
Eq.|8| with the actual ones. 

There are no GEANT interfaces to the production models other than FLUKA. 
Different production models can be studied by re-weighting the neutrinos according 
to the production cross section of its parent. Such a weighting procedure, developed 
by Mark Messier [IT]], exists only for the parent 7r + . 

The correlation matrix M was re-computed for the neutrinos produced by the 
pions only and leaving aside the neutrinos produced from kaon and muon decays as 
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well as neutrinos produced in proton-induced showers in the final absorber. This 
matrix was subsequently applied to the near detector spectrum, re-weighted for 
different models and the results were compared with the far detector spectrum in 
these models. Results are shown in Fig.|18[ 
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Figure 18: Upper figure: ratios of far detector spectra predicted by various pro- 
duction models to those predicted from the corresponding near spectra using the 
correlation matrix M calculated with the GFLUKA model in ow energy beam. 
Lower figure: a comparison of the Far/Near ratios in different prodution models 
over the nominal GFLUKA Far/Near ratio. 

Far detector spectrum is predicted with the accuracy of the order of 2% for 
the component of the neutrino beam produced by the focused part of the hadron 
beam. The prediction is accurate to ~ 5% for neutrinos energies up to ~ 12 GeV, the 
sensitivity of the method to the production spectra somewhat reduced in comparison 
with the 'double ratio' method. 
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4.3 Production Model Dependence: Medium Energy Beam 

The procedure described in the sectionIO was repeated for the case of the medium 



energy beam. The results are shown in Fig. 19 
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Figure 19: Upper figure: ratios of far detector spectra predicted by various pro- 
duction models to those predicted from the corresponding near spectra using the 
correlation matrix M calculated with the GFLUKA model. Medium energy beam 
option. Lower figure: a comparison of the Far/Near ratios in different prodution 
models over the nominal GFLUKA Far/Near ratio. 



As before, the prediction is good to ~ 2% for the entire region of peak intensity 
and it remains good to ~ 5% for energies up to ~ 12 GeV . The sensitivity to the 
production model details is, again, reduced in comparison with the 'double ratio' 
method. 
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4.4 Sensitivity to the Details of Modeling of the Neutrino 
Beam 



Functions F n (ri, 6,p) and P{r^ 9,p, z) used to compute the correlation matrix matrix 
M depend on the momentum and spatial distribution of hadron in the decay volume, 
hence they depend on the geometry and focusing characteristics of the beam forming 
elements. How sensitive is the resulting matrix to the details of the understanding 
of the beam line? 

Two programs are currently in use to simulate the neutrino fluxes: 

• GNUMI is a GEANT-based full simulation of the production and focusing of 
particles. It includes scattering and re-interactions of the produced hadrons. 

• PBEAM is a much faster, albeit somewhat simplified, intended for optimiza- 
tion of the beam line design. 

We have used these two programs corresponding to somewhat different beamline 
configurations: geometry of the target cave, details of the magnetic field in the horn, 
thick window of the decay volume are among the known differences between them. 
It is not surprising, therefore, that the shape of the neutrino spectrum resulting 
from these two programs are not identical, as shown in Fig.[20| In a similar fashion 
these two programs yield different spectra for the far detector. 

From the spectrum at the near detector the far detector spectrum can be calcu- 
lated using the Eq.[|. A comparison of this prediction with the actual far detector, 
PBEAM-calculated, spectrum is shown in Fig. [21] together with the similar compar- 
ison for the 'double ratio' method. Despite the major differences in the simulation 
programs and the resulting difference in the predicted neutrino spectra between the 
program used to derive the correlation matrix and the program used to 'analyze' 
the data the prediction is good to a level ~ 5% or better. Again, the 'double ratio' 
is more sensitive to the differences between the program than the matrix method. 

4.5 Distribution of Pions in the Decay Volume 

As shown in previous sections, the difference between the near and far detector 
spectra is directly related to the finite decay pipe size and to the distribution of 
decay vertices inside the decay volume. 

A robustness of the prediction of the far spectrum can be investigated by arti- 
ficially enforcing major deviations from the spatial distributions produced by the 
beam simulation programs. While such a procedure does not correspond, probably, 
to any realistic imperfections of the beam line, it nevertheless offers an insight into 
a sensitivity of the method. 

A radial beam profile can be modified by applying a weight factor of the form 




(10) 
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Figure 20: Up: near detector spectra as predicted by PBEAM and GNUMI programs 
Down: a ratio of the near detector spectra in these two programs 
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Figure 21: Up: ratio of the far detector spectrum derived using GNUMI-based 
correlation matrix to the actual one as predicted by PBEAM. Down: the same ratio 
for the far detector spectrum predicted using a 'double ratio' method. 



where r = 50 cm corresponds to a half-radius of the decay pipe and a = 0.2, 0.5, —0.2, 
is the magnitude of the distortion. 

Application of the standard transformation to the near detector neutrino spec- 
trum corresponding to the shrunk down or blown-up hadron beam leads to a pre- 
diction of the far detector spectrum with the accuracy shown in Fig .|2"2|. 

Prediction of the shape of the neutrino spectrum in the far detector is not very 
sensitive to the details of the transverse beam profile. Shrinking and expanding the 
transverse profile of the hadron beam by ±20% affects the predicted far detector 
spectrum by 2% or less. Even more drastic change of the profile by ±50% leads to 
a modest < 5% change in the predicted spectrum. 

Variation of the distribution of the decay vertices along the beam axis is poten- 
tially more serious source of the systematic error, given the proximity of the near 
detector. A sensitivity to the knowledge of this distribution can be evaluated in a 
similar manner by re-weighting the decay points with 

w z (z) = 1 + a ■ [z - z \/zq (11) 
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Figure 22: Up: ratios of far spectra for different beam profiles to those predicted 
from the corresponding near spectra using the same "nominal" M matrix. Beam 
profiles were defined by changing parameter a in formula (|T0l) . Down: an equivalent 
comparison of Far/Near ratios over the nominal Far/Near ratio. 



where zq = 300 m and a is variable parameter, as before. 

Results, shown in Fig.|23"|, indicate that the knowledge of the longitudinal distri- 
bution of the decay points to ±20% is needed to keep the systematic error within 
the bounds of the statistical accuracy. It is not clear what effects in the beam line 
may cause such a significant re-shuffling of the decay points. One possible source of 
such an effect would be an additional component of the hadron beam consisting of 
well collimated pions thus traveling much longer distance along the decay volume. 
Such a component can be created, for example, by a proton beam scraping inside 
of the collimators upstream of the target. High energy pions created in such inter- 
actions and focused by magnetic horn would have small angular divergence (owing 
to the large longitudinal momentum component) and hence they would produce 
disproportionally large neutrino flux in the near detector. 

To keep the systematic errors under control it is therefore very important to 
minimize the amount of materials which can be intercepted by the proton beam 
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Figure 23: Up: ratios of far spectra for different pion decay vertex distributions 
to those predicted from the corresponding near spectra using the same "nominal" 
M matrix. Vertex distributions were modified by changing parameter a in formula 
flnp. Down: an equivalent comparison of Far/Near ratios over the nominal Far/Near 
ratio. 



upstream of the target. It is perhaps possible to constrain the actual distribution of 
the hadron in the decay pipe , and thus reduce the systematic error, by monitoring 
the energy deposition along the length of the decay volume. 

4.6 Horn Current 

The event rates and their energy distribution in the far and near detectors, as well 
as the far/near ratio depends on the strength of the focusing elements. These effect 
was investigated using samples of events generated with GNuMI the horns current 
changed from the nominal value of 200 kA to the values of 180 kA and 220 kA. The 
near detector spectra were used to compute the far detector spectra using the matrix 
M corresponding to the nominal horn current. 

Fig.|24] shows that to maintain the neutrino flux prediction with the accuracy 
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below 5 % it is necessary to monitor the horn current to 5% as well. In addition 
to the monitoring of the peak current for the horn, the timing of the current pulse 
must be monitored as well. In a particular case of the NuMI horn with 1.7 msec 
pulse it is required that the timing of the pulse is monitored to better than 150 fisec 
to maintain the required level of the systematic error. 
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Figure 24: Up: ratios of far spectra for different horn currents (red: I = 0.18 MA, 
green: I = 0.22 MA) and those predicted from the corresponding near spectra using 
the same "nominal" M matrix. Down: an equivalent comparison of Far/Near ratios 
over the nominal Far/Near ratio. 



4.7 Horn Displacements 

Position of the magnetic horns defines edges of the beam acceptance and the position 
(in the phase space of the produced particles) of the boundary between the focused 
beam and the bare target beam. Small displacement of the horn causes some bin of 
high momentum pions to be focused and increase significantly their contribution to 
the neutrino flux, at the expense of the small fraction of neutrinos at slightly lower 
energy. Displacement of the first horn will produce also change of the phase space of 
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particles focused by the second horn, too. The resulting change of the event spectra 
observed in the near and far detectors is shown in Fig.p5. 
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Figure 25: Ratios of the event spectra in the near (top) and far (bottom) detectors 
observed with the horn 1 displaced by 2 mm to the spectra with horn 1 in the 
nominal position. 



Due to the finite radial extent of the decay volume the additional component of 
the neutrino flux shows up at slightly different energies in the near and far detectors. 
As a result, the sensitivity of the 'double ratio' method to the horn displacement is 
somewhat enhanced. 

The far detector spectrum predicted with the help of the correlation matrix M, 



Eq|| is much less sensitive to the horn displacements, as illustrated in Fig. 26 



5 Further Improvements: Large Near Detector 

The correlation matrix M defined in Eq.fi] relates the near and far detector spectra 
averaged over the detectors volume. The far detector, in all practical cases, is located 
so far from the neutrino source that there is no variation of the neutrino spectrum 
over the area of the detector. 
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Figure 26: Ratios of the predicted and observed spectra in the far detector when 
horn 1 is displaced laterally by 2 mm. Top figure is for the prediction using the 
matrix M corresponding to the nominal horn position, bottom figure is for the 
double ratio' method. 



The transverse size of the MINOS near detector is small compared to the typical 
beam size of the low energy beam, hence the near detector spectrum is practically 
identical over the entire fiducial volume. 

In case of a near detector large enough to detect spatial variation of the neutrino 
flux the correlation matrix can be further refined to take advantage of the additional 
information. This might be the case, for example, for the K2K experiment. 

A correlation matrix M^a can be calculated by carrying out the integral ^ us- 
ing the weighting function WdA( z , r dec, Q,P, z n, E n ) describing the probability that a 
pion(kaon) with momentum p and the angle 9 decaying at the position z along the 
decay volume at the radius r dec will produce neutrino with energy E near inside the 
detector area dA in the near detector. 

The far detector spectrum can be computed by integration over the near detector 
area of the predicted contributions, as implied by the spectrum of events detected 
in the area elements dA: 
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(F U F 2 ,--- ,F n ) = j 
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dA 



(12) 



where the vector N d A = (Ni,N 2 , ■ ■ ■ ,N n ) dA is the spectrum of the events de- 
tected in the area dA of the near detector. 



6 Conclusions 

A prediction of the neutrino flux in the far detector from the one observed in the near 
detector is very robust. This prediction is primarily determined by the geometry of 
the beam line and the properties of the focusing elements. 

The overall event rate can be predicted with an accuracy better than 2%. Sys- 
tematic error on the shape of the energy distribution of neutrinos are of the order 
of 2 — 3% for the main part of the flux (corresponding to the beam created by the 
well focused fraction of the hadron beam) and of the order 5 — 10% in the tails. 

The sensitivity of the far detector flux prediction is maximal in the energy regions 
corresponding to the edges of the acceptance of the magnetic horns. In a particular 
case of the low energy NuMI beam it is a region around 5GeV and 8 — lOGeV, 
corresponding to the acceptance edges of the horn 1 and horn 2, respectively. 

An improved prediction method, utilizing the correlation of the energy spectra, 
reduces the sensitivity of the far detector spectrum prediction, especially to the ef- 
fects producing localized distortions of the spectrum. This method enables a reliable 
prediction of a neutrino flux at different distant locations, based on the same near 
detector measurement. 
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